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ABSTRACT 

Following the discovery of a black hole (BH) with a mass of 10^-10^ Mq in a 
starburst galaxy M82, we study formation of such a BH via successive merging of 
stellar-mass BHs within a star cluster. The merging has a runaway characteristic. 
This is because massive BHs sink into the cluster core and have a high number 
density, and because the merging probability is higher for more massive BHs. 
We use the Smoluchowski equation to study analytically the evolution of the BH 
mass distribution. Under favorable conditions, which are expected for some star 
clusters in starburst galaxies, the timescale of the runaway merging is at most of 
order 10^ yr. This is short enough to account for the presence of a BH heavier 
than 10^ Mq in an ongoing starburst region. 

Subject headings: black hole physics — galaxies: starburst — galaxies: star 
clusters 



1. INTRODUCTION 

While black holes (BHs) arising from stellar evolution have masses of orders 10°-10^ Mq, 
those found in galaxy nuclei have masses of orders 10^-10^ Mq (Kormendy & Richstone 1995; 
Kormendy & Gebhardt 2001). The possible missing link, i.e., an intermediate-mass BH, was 
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recently discovered in a starburst galaxy M82 as a source of compact X-ray emission (Kaaret 
et al. 2001; Matsumoto et al. 2001; see also earlier references therein). The observed strong 
variability implies that the source is an accreting BH. The observed luminosity of 10^^ ergs s~^ 
implies that the BH mass is greater than 10^ Mq if the emission is isotropic and its luminosity 
is below the Eddington limit. The BH lies at the 2 yum secondary peak (Matsushita et al. 
2000), which is an unusually active site of ongoing star formation. Hence the BH was formed 
during the ongoing starburst and is less massive than a so-called super star cluster, i.e., a 
building block of a starburst region. The typical mass of a super star cluster is of order 
IO^Mq. 

The formation of an intermediate-mass BH is attributable to successive merging of 
stellar-mass BHs within a super star cluster (Taniguchi et al. 2000; Matsushita et al. 2000; 
Ebisuzaki et al. 2001; Miller & Hamilton 2001). The cluster sinks into the galaxy center due 
to dynamical friction and upon evaporation releases the intermediate- mass BH. If there are 
more than one such BHs, they merge with each other. Thus an intermediate-mass BH could 
eventually evolve to a supcrmassivc BH and is crucial to studying a possible evolutionary 
connection between a starburst and an active galactic nuclei (Mouri & Taniguchi 2002 and 
references therein). The latter harbors a supermassive BH as the central engine. 

Since the 2 yum secondary peak of M82 has a starburst age of order 10^ yr (Satyapal 
et al. 1997), the formation of the intermediate-mass BH should have occurred within this 
short duration. There is a possibility of runaway merging, which is known for collisions of 
BHs and those of normal stars (Quinlan & Shapiro 1989; Lcc 1993; Portcgics Zwart et al. 
1999). Massive objects segregate into the cluster core due to dynamical friction, increasing 
the merging probability. The probability itself is higher for more massive objects. Their 
merging products are even more massive. Thus the growth rate is higher than that expected 
from the merging rate at the initial time. 

However, runaway merging is a nonlinear process and hence has been studied only for 
specific cases using numerical methods, e.g., iV-body simulation. Theoretical understanding 
of the underlying physics is insufficient. Also, the existing numerical results are not applicable 
to formation of an intermediate-mass BH. A more general analytical approach, albeit less 
realistic, is desirable. 

We make an analytical analysis of the runaway merging. In our idealized model, BHs of 
a single mass mo exist initially with a number density Uq and a velocity dispersion vq. These 
BHs are assumed to merge with each other via two-body interactions, i.e., energy loss due 
to gravitational radiation (§2). The BH mass distribution is studied using the Smoluchowski 
equation, a master equation describing the evolution of number densities of various-mass 
particles that merge with each other via two-body interactions. We obtain an analytical 
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constraint on the timescale of runaway merging (§3), and apply the result to BHs in a super 
star cluster with mo = 30 Mq, no — lO^pc"^, and Vq — lkms~^ (§4). 



2. MERGING VIA GRAVITATIONAL RADIATION 

When two BHs become close, energy loss due to gravitational radiation could exceed 
the orbital kinetic energy, and a binary could form. The binary immediately merges through 
subsequent gravitational radiation. For this process, we estimate the cross section and then 
the rate coefficient. 

We consider an interaction of two BHs where the BH masses are rrii and rrij, the initial 
relative velocity is v^eh and the impact parameter is b. If the two BHs pass away narrowly, 
the orbit is close to parabohc. We thereby assume that the eccentricity e is shghtly greater 
than 1, and consider the leading term for e — > 1. This is equivalent to an assumption of 
significant gravitational focusing. The distance of closest approach rmin is 

' TTiin — 



2G{mi + rrij) 

For an unperturbed parabolic (e = 1) or hyperbolic (e > 1) orbit, the quadrupole formalism 
gives the total energy of the gravitational radiation 5E aX, Newtonian order as 

^E = T^-. ^2 9{e), (2) 

min 

with g{e) = 4257r/(32\/2) at e = 1 (Turner 1977). The condition for the two BHs to pass 
away is 

2mi + mj '""^ ^ ' 

Equations (l)-(3) yield the minimum impact parameter femin for the BHs to pass away. Then 
we obtain the merging cross section (Tmer = '^^^m ^ 

"-"'Hi^j ■ 



Quinlan & Shapiro (1989) derived the same formula with a somewhat different reasoning. 
The cross section cTmer is larger than that corresponding to the Schwarzschild radius (Tsch — 
7r(2Gm/c^)^ ~ ('yrei/c)^^^^Cmer- Thus our Newtonian approximation is sufficiently accurate. 
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Since the merging cross section (Jmer is smaller than that for gravitational scattering, 
(7sca — 7r(Gm/i'^gi)^ ~ (c/vreO^^^^Cmer, the BH motions are assumed to be in thermal equilib- 
rium as a zeroth-order approximation. The equipartition of the kinetic energy is assumed to 
be achieved as ^ ^ 

-rriivf ^ -itioVq for any rrii. (5) 

Here Vi and Vq are three-dimensional velocity dispersions for BH masses rrii and mo, respec- 
tively. The probability distribution of the initial relative velocity Vrei is given by 

^/ X / x\ . , 3/;/,/;/ , (',•^,,1 

p{x) = ( — ) exp -- with x = - „ . 6 

Equations (4) and (6) yield the average (vreiCmer) as the rate coefficient R^^j for merging of 
BHs of masses and mf 

AC^ml fVo\ -11/7 / rrnmj \ ^^^^^ f mi + rrij \ ^^^^ 

with A — r(5/7) (27r)ii/i^ 3^/^ 85^/^ ~ 33.33. However, in this estimate, mass segregation 
has not been taken into account. Massive BHs sink into the core of the star cluster and have 
a high number density even if their total number is small. The scale height of BHs of a mass 
rrii is proportional to (multimass King's model; Sigurdsson & Phinney 1995). Thus 

mass segregation enhances the merging rate by a factor {rriimj /m^Y^'^ . The rate coefficient 
is then proportional to {mirrij/miy^^'^ instead of {mirrij/miy^^^^. 

The expression (7) is too complicated to be studied analytically using the Smoluchowski 
equation. Wc use the inequality rrii + iT^j ^ lirniVrijY/'^ and obtain a simplified expression 
for the merging rate coefficient: 

with B = r(5/7)7rii/i^2i°/^3i/2 852/7 ~ 52.04. The exponent A is 39/14 or 81/14, respec- 
tively, if we ignore or incorporate the mass segregation. 

We consider two-body interactions alone and thus ignore three-body interactions. They 
result in BH binaries. The BH binary undergoes subsequent interactions with other BHs and 
increases its binding energy. At each of the interactions, a fraction of the change in binding 
energy is converted into the translational kinetic energy. The BH binary eventually merges 
by emitting gravitational radiation or escapes from the star cluster (Sigurdsson & Hernquist 
1993; see also Portegies Zwart & McMillan 2000). Although a numerical calculation based 
on a reliable BH mass function is required to draw a definitive conclusion, we expect that 
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the former occurs much more often than the latter in a massive cluster with a large escape 
velocity if the BHs are massive enough 10 M©; Miller & Hamilton 2001). The merging 
timescale is short for a binary of massive BHs. The escape of a massive BH binary is 
due to a collision with a massive BH. Before such a rare collision, the BH binary would 
merge. Our present model would accordingly underestimate the merging rate. This possible 
underestimation is not serious because we are to obtain an upper limit on the timescale of 
the runaway merging. 



3. SMOLUCHOWSKI EQUATION 

Suppose that particles of a single mass rriQ exist with a number density Hq at the initial 
time t = 0. They merge with each other via two-body interactions. The merging of particles 
of masses nii and nij has a rate coefficient Rij and results in a particle of a mass + mj. 
We consequently have particles of masses = moi {i — 1, 2, 3, ...) and number densities rii. 
If nondimensional terms are defined as 

■hi^—, Rij^^^, and i ^ noRo,ot, (9) 
1^0 -Ko,o 

the evolution of the number densities is described by the Smoluchowski equation: 

-^^\Y1 ^hkfijfik -rii^ Rijhj. (10) 

j+k=i 3=1 

Here the first sum accounts for the increase of hi due to merging of particles satisfying the 
condition mj+mk = rrii while the second sum accounts for the decrease of hi due to merging 
of particles of a mass rrii with those of any mass. The initial condition is ni = 1 and = 
for i > 2 at f = 0. 

The Smoluchowski equation is based on a mean-field approximation and thereby ignores 
spatial fiuctuations, e.g., mass segregation. It is nevertheless possible to incorporate conse- 
quences of the fiuctuation into the rate coefficient Rij. The number densities are then 
regarded as those averaged over the fluctuation scale for particles of a mass mo. 

The analytical solution of equation (10) is known only for Rij = 1, {i+j)/2, and ij (see 
Hayashi & Nakagawa 1975). We study the solution for Rij = {ij)^^"^ with A > 2 and relate 
it to the solution for i?j j > (ij)^''^, which corresponds to our rate coefficient (8), using the 
^-th moment: 

oo 

M, = Y,^%. (11) 

i=l 
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The total mass density and average mass of the particles are given respectively by Mi and 
M2/M1 (Hayakawa & Hayakawa 1988). Prom equation (10), the moment evolves as 

^ = ij]n.nA.[(^ + Jr-^^-/]- (12) 

We have dMi/dt = 0. which means the conservation of the total mass, Mi = 1. Multiplying 
both sides of equation (12) for Rij — (ij)^^'^ and £ = 2 by M2~^ gives 

— < -(A - 1) for A > 2. (13) 

dt 

Here we have used the mass conservation Mi = 1 and also a relation originating in the 
Schwarz inequality: 

^ >(%i^) for A > 2. (14) 

The equahty holds at A = 2. Equation (13) has the solution: 

1 



M2 > 



l-(A-l)t 



for A > 2. (15) 



Thus the average mass M2/M1 — M2 becomes infinity at a finite time frm < (A — 1)~^, owing 
to emergence of infinite-mass particles. We regard as the timescale of runaway merging, 
which is shorter than the timescale t = 1 corresponding to the initial merging rate UqRq^q 
(eq. [9]). Even if Rij = {ij)^^'^ is replaced with Rij > (ij)^^^, the solution (15) remains the 
same. 

The total number density of the particles is given by the zeroth moment Mq. Hence it 
would appear natural to adopt Mi/ Mq = Mq'^ as the average mass. If both sides of equation 
(12) for Rij > (ij)^/'^ with A > 2 and £ = are multiplied by Mo^"^ it follows that Mq 
becomes zero at a finite time ^ 2(A — 1)~^ However, before the time moment 
equation Mi = 1 becomes invahd (Hendriks, Ernst, & Ziff 1983). This is the reason why we 
have adopted M2/M1 as the average mass. 

The emergence of infinite-mass particles occurs in more general cases. If the rate coef- 
ficient Rij scales as Rai,aj = oj^Rij = Oj^Rj,ii infinite- mass particles emerge within a finite 
timescale for A > 1 (Hayashi & Nakagawa 1975). Lee (1993) used this property to interpret 
runaway merging of BHs found in A'"-body simulations. Since the timescale of the runaway 
merging had not been estimated analytically, we have obtained the analytical constraint (see 
Hendriks et al. 1983 for similar discussion in terms of sol-gel phase transition). If the rate 
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coefRcient behaves as Rij ~ i'^j" with > 1 for j — > oo, infinite-mass particles emerge and 
high-order moments become infinity just after the initial time t = (van Dongcn 1987). 
This is the case in our rate coefficient (8). We have ignored moments with £ > 2 because 
the Smoluchowski equation (10) is merely used as an idealized model in the present work. 



4. Discussion 

The runaway merging of BHs is caused by two mechanisms, which are equally important. 
First, the merging rate Ri^j is intrinsically large for massive BHs (eq. [7]). This is because 
the energy loss due to gravitational radiation 6E is large. Massive BHs have low velocities 
(eq. [5]). If the initial relative velocity fj-ci is low or the BH masses and rrij are large, the 
gravitational focusing is significant and the distance of closest approach rmin is small (eq. 
[1]). If the distance rmm is small or the BH masses rrii and rrij are large, the energy loss 6E 
is large (eq. [2]). Second, low- velocity massive BHs sink into the cluster core and have an 
enhanced number density. The first and second mechanisms correspond to A = 39/14 and 
3, respectively. They are comparable and individually exceed the critical value for the onset 
of runaway merging, A = 1 (§3). 

Equations (8), (9), and (15) yield a constraint on the timescale trm of the runaway 
merging as 

trra < T T ^r^. o W 



1 BC^mlrio \ c J 

mo \ ( TiQ \ ( Vq \ J 4 X lU yr no mass segregation 

X 



30 M©/ \10^pc \lkms ^/ \ 2 x 10''' yr mass segregation. 

The mass mo = 30 M©, the number density Uq = lO^pc"^, and the velocity dispersion 
Vq = lkms~^ provide the most favorable conditions for the runaway merging among those 
expected for BHs in actual super star clusters (see below). We accordingly obtain an upper 
limit of order 10 ''yr on the timescale trm- This upper limit is comparable to the starburst age 
observed for the 2 //m secondary peak of the starburst galaxy M82 (§1). Since a starburst 
duration is of order 10^ yr (see de Grijs, O'Connell, & Gallagher 2001 for the case of a fossil 
starburst M82 B), it is generally possible that an intermediate-mass BH is formed in an 
ongoing starburst region. 

The BH mass mo = 30 M© has been adopted because stars heavier than 30-40 
evolve directly to BHs with no supernova explosion (Brown & Bethe 1994; Fryer 1999). The 
BH mass is equal to that of the progenitor. Although the less massive stars could evolve 
to BHs, they explode as supernovae and blow off most of their mass. A low-mass BH is 
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easily ejected from the star cluster by three-body interactions and does not contribute to 
the intermediate-mass BH (§2). 

The number density no — lO^pc"^ and the velocity dispersion vq — 1 kms~^ of the BHs 
have been adopted in the light of a nearby super star cluster, R136 in the Large Magellanic 
Cloud (see Portegies Zwart et al. 1999). The central mass density of R136 is lO^Mopc"^. 

Wc assume the Salpeter initial mass function dn/dm oc m~^'^^ with the cutoffs at 0.1 and 
100 Mq, which yields the number fraction 0.04% for stars heavier than 30 Mq and the average 
mass 0.4 Mq. Wc also assume the multimass King's model (§2). Then the above mass density 
corresponds to the number density no ^ 10^ pc^"^ for BHs heavier than 30 M0. Here wind 
mass loss from stars has been ignored as a zeroth-order approximation. The mass loss rate 
observed for a 30 Mq star is 10~^ Mq yr~^ (Becker 1999). Since the lifetime of a 30 Mq star is 
10^ yr, the total mass loss is at most a few times Mq. On the other hand, the stellar velocity 
dispersion of R136 is estimated to be about lOkms"^, which corresponds to the velocity 
dispersion vq ~ Ikms^^ for BHs of a mass 30 Mq. If runaway merging occurs among stars 
before they evolve to BHs (§1), the number density uq could be higher and hence the upper 
limit on the timescale trm could be lower. 

The timescale t^is of mass segregation is less than that of the runaway merging. If 

the average stellar mass m is 0.4 Mq, the mass density p is IO^MqPC"^, and the velocity 
dispersion v is lOkms""^, the star cluster has the relaxation timescale t^i^ — / (irG'^mp) = 
4 X 10^ yr, which yields tms — (^/?7T.o)^rix = 5 x 10^ yr for BHs of a mass mo = 30 Mq (Spitzer 
1969). 

Super star clusters have ranges of structure parameters. Existing few observations of 
super star clusters other than R136 tend to imply vq > lkms~^ (Ho & Filippcnko 1996; 
Smith & Gallagher 2001). Observations of globular clusters, i.e., likely descendants of super 
star clusters, tend to imply Uq < 10^ pc^^ (Pryor & Meylan 1993). In these cases, the 
timescale trm of the runaway merging is long (eq. [16]). The formation of an intermediate- 
mass BH in an ongoing starburst region occurs only under rare favorable conditions, i.e., 
high number density Uq and low velocity dispersion vq of stellar-mass BHs. Also, the star 
cluster as a whole has to be sufficiently massive in order to provide a sufficient number of 
the stellar-mass BHs. In fact, only one intermediate-mass BH heavier than 10^ Mq has been 
found among more than 100 super star clusters of the starburst galaxy M82, although there 
might exist some intermediate-mass BHs that are not accreting the gas and hence are not 
observable. 
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